“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1971-09 


The formation of hydraulic jumps ina 
rotating, continuously stratified system. 


Kurth, Robert Paul 


Monterey, California ; Naval Postgraduate School 
http://ndl.handle.net/10945/15726 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


i 2 ae 


THE FORMATION OF HYDRAULIC JUMPS IN A 
ROTATING, CONTINUOUSLY STRATIFIED SYSTEM 





Robert Paul Kurth 











by 


+ 


SE ETE, CB CS ah SE 


Roberti Paul Kurth 


Re IT. Wilthams 


Co=-Thesis Advisors: 
Gade blattines: 








THE FPORMATION/SOF Hy DRAUL C2 UiiEs 
IN TAL ROTA TING. 
CONTINUOUS y SiRA err peewee! ei 

Pepreciaiet h(t 


cS eaniiasiaiaaaanel ce mee Do ee ee ee eee 








meee 
NE PE Se A CD pA EE EES ARE OS TR Ne i ee ee ee 


Approved for public aclease; distubution wiloncted. 








ihe Formation of Hydraulic dumps anvashotacing- 


Continuously Stratified System 


by 


Robert PawliGiees: 
Lieutenant, United States Navy 
B.S.E., Kansas State Teachers College, 1966 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN METEOROLOGY 


from the 


Naw AL POs tT CRADUATE SCGHeoor 
September 1971 








ABSTRACT 


The formation of hydraulic jumps ina rotating, continuously 
stratified, Boussinesq fluid is investigated. A scale analysis indicates 
that jump-formation cases should be separated from non-jump cases 


baethe curve Hl = ARo* 


where Ro and Fare the initial Rossby and 
mooude mumibers, respectively. Numerical solutions verify the curve 
up to F = 0.3 and show that A is inthe range 6.0 to 7.5. The solutions 


for F = 1.0 are examined in considerable detail. The results of this 


paper extend the work of Houghton and Wiiliams and Hori. 
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lL, ANTROBUCTION 


The formation of hydraulic jumps has been shown for a one-layer, 
homogeneous fluid for both non-rotating and rotating cases. This 
eeperiment willl consider a one-layer, continuously stratified fluid by 
first investigating jumps in a non-rotating coordinate system and then 
extending the study to the rotating system. 

In a non-rotating one-layer fluid a hydraulic jump with a discon- 
tinuity in the velocity can arise from an initial state having any scale. 
If the system rotates, there may bea limiting initial scale beyond 
which jumps cannot form as shown by Houghton [1969] and Williams 
eecemon: | ORC}. Aneta: pemtohbavicr in exmmeqgred in a continucwSiy 
Stratiiied atmosphere, In non-dimensional terms this corresponds 
toa critical Rossby number (Ro - U/fL, where U is the initial velocity 
scale, Lis the initial space scale and fis the Coriolis parameter) 
below which jumps are excluded. This process also depends on the 
roude number (F = U/[D(20%,!40/a2) We icre D is the vertical scale 
and 6, is the mean potential temperature). In the F-Ro plane there 
should exist a critical curve that would separate the jump cases from 
the non-jump cases as has been shown for a homogeneous fluid by 
Houghton [1969] and Williams and Hori [1970]. In this experiment it 
was shown that a critical curve exists and is valid for time periods of 


e 


= 
order L/[D(g6. d@,/dz) ]. For longer time periods, however, jumps 








may occur in the non-jump region, indicating that the critical 'curve' 
is actually a zone where the jump time increases rapidly. 

There are several possible applications of this work in meteorology 
and oceanography. Rotation does have an effect on mesoscale motions 
in the atmosphere, but the influence is not dominant as it is in large= 
scale motions. Tepper [1952] has proposed that squall lines are 
modified hydraulic jumps. Freeman [1948], Abdullah [1966], and 
Houghton and Kasahara [1968] have proposed other mesoscale hydraulic 
analogies. Sasaki [1959] and Ogura and Charney [1962] have studied 
squall lines in continuously stratified atmospheres. 

Internal waves in the ocean are affected by rotation, and the break- 
ing of these waves should be described by this model. Rapidly moving 
frontlike temperature disturbances have been observed in the ocean 
[Gowans, 1969] and these may be internal waves modified by rotation. 

In this experiment the Boussinesq equations are given and the 
modeling relations are introduced inl]. In the model the time-dependent 
quantities are functions of x and z only. The quantity f is assumed to 
be constant for each independent case study. The finite difference 
approximations are presented in Ill. A scale analysis is presented 
+n IV which is useful in the discussion of the numerical solutions. ive 
initial conditions are given in V and the numerical solutions are pres- 
ented and discussed in VI. The initial state consists ofa sinusoidal 
disturbance (amplitude dependent on F} which is superimposed on the 


basic state. After a sufficient period of time, the sine wave 1s distorted 








in sucha way that a hydraulic jump is formed. In VII conclusions and 


recommendations for further investigation are given. 








II. BASIC EQUATIONS 


In this experiment the Boussinesq equations are employed and 
the domain is bounded by two rigid horizontal planes. The compress- 
ibility of the atmosphere is neglected in the Bossinesq approximation. 
This should not be of qualitative importance since the density scale 
height of the atmosphere is much larger than the vertical scale of the 
typical jump zone. The hydrostatic approximation is also used in 
this experiment and it should be good as long as the horizontal scale 
is much larger than the vertical scale [Waa LOA Se paar one 

The Boussinesq equations which were obtained by Ogura and 
nes 2 las 


Soret. Oba ag [aeons £ yon 07 he srvitte rae a7 : mies, Bee Fa * —,<- 
a ahiaslpo L Jah re VAT 4 fo gr al Vy Lac LPL We APS @eh WresaUad Caosaw ee ae oe oe 


rotation is included: 


aV en, Ooty, ev, 
<= + 0(VV) + 5g (wV) +VP+Fkx V0, (II. 1) 


rt Wee a (w6) =O, (1. 4) 


aa.) 
VV tyz =O ana aie) 


3P g © (II. 4) 











ie constant reference potential temperature, 
Po = constant reference pressure, 
= Rep C _, 
K i 
8 = T(P,/P). - 6, = departure of the potential temperature 


aces [shan einayal 

6 = C,0,(P/P,) + gz. 

Heating and friction have been neglected in these equations even though 
they may be important under certain conditions. 

Phipguahitties are independent Of @y in this model. Tne Goro. 
parameter, f, is assumed to be constant throughout the model but is 
varied for individual cases. A Cartesian coordinate system is intro- 
duced with x increasing toward the east and y increasing toward the 
north. The x and y velocity components are u and v respectively. 


For this model, Equations (II. 1)-(II.3) can be written in the 


following form: 





au (uu) d(uw) 39 = 2 (ii 
ae Tt Te re ae 2 ae fy =O) 


dv y(uv) . d(vw) _ 
ea Pag Pe =O, 


dO y(u 6) »(w 8) - 


oe See eee (II. 
3 T Yx 2 dZ = 0 ano 


=~) 
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The boundary conditions are 

Ww = O. 2 520.) (TL293 
where D is the distance between the horizontal plates. 

The integral of Equation (II. 4) gives the departure of $ from its 
vertical average, and this average can be obtained as follows. The 


following averages are defined: 


= wl 0 
( errai ihe dz, ¢ dep )dz, (II. 10) 


caeere 2 ill Ibis the periodicity of the fields inx. If Equation (1176) is 


averaged vertically and (II. 9) employed 


a Glens = Cae (ts) 





Eumagjarly, the vertical average of (11.5) is 


re ase. ST ae (II. 12) 


If this equation is in turn averaged with x and if periodicity in x is 


imposed, Equation (1.12) can be written 


d((%)) 


at 





~ fa.7)y = O, (II. 13) 


The quantity (9) represents the net mass Mux in the y direction. Wi 


((¥)) were different from zero, (=) could not be zero if a wall 


were present at some value of y. Since this would violate the boundary 


conditions [Equation (II. 9)], the following is a requirement in this 


See miner: 


Lal 








((7)) = 0. (II. 14) 


itthen follows directly from (II. 13) that 


J((u)) 
dt 





= Q, (II. 15) 


dud de) 


Now from (Il. 11) which together with (II. 15) 


st = =k 
reduces (II. 12) to d¢¢d) - om MCuuy 4 Sev). GL, te) 
2X aX 


The corresponding equation is obtained by averaging Equation (iI.6) in 
x and z and by introducing (IJ. 14), i.e., 

f<cws =0. (11. 17) 
Equation (II. 17) shows that the x- and z-averaged zonal wind is zero. 
ge o.4), {12.15) ae 


dite MY US LL GTS ‘sc. oy A (= 


(II.6) can be rewritten in the forms 


bu + s(uUu) a “Oe 2 (uw) ative -(v)) 


a 





: _ dx oz (II. 18) 
4: (9 -¢#))=O and 

Vv dCuv) d(wv) 

a ESE ari le eS (II. 19) 

ye OUT a a8 Su = 0, 


By integrating the hydrostatic equation [Equation (Il. 4)] the following 


expression for O- <4) is obtained: 


9 -<¢) = 3 ['s dz = (fe x) |. (II. 20) 


The following energy integral is an invariant for this model: 


pe 1d uy) + € »)| 2 : 3((22)). (II, 21) 


I 











iit, FINITE DIFFERENCE APPROXIMA EONS 


Equations (II.7), (11.8), (II. 18) and (IJ. 19) were solved numerically 
Pyeantroducing finite differences in x, zandt. The solutions were 
obtained by a simple marching process, but if the horizontal variations 
were two-dimensional, Poisson's equation would have to be solved at 
each time step [Smagorinski, 1958; Hinkelmann, 1959]. 

The vertical domain was divided into 2K layers of equal thickness, 
and the surfaces separating the layers were denoted by k. The arrange- 
ment of variables on the bounding surfaces is shown in Figure l. The 
eiverpence oF tic vertical flix ofitme Scalar at 2a Z is approximated 


ads 


E a : Wk +t (Ss Whe! (Si + Seay) ait. 4) 
k 


dz zZ (Zia ay) 





k= Ars ave © 
f= K-| wa 
rd U7 
poe VECO. Ieee 
2 
Dr ee ee Dc 
kes ee | 
keO w= 0 Z=0 


Fig. 1, Ayrangement of variables im Z. 


Ws 








The finite difference approximation to the integral in Equation (II.20) 


used is 


JF 


{e = [G2 (2,,-za) $C, k>s. (III. 2) 


ate 


The value of C fret 1S immaterial since it will disappear when the 
vertical average is removed. The finite difference integral of Equation 


(If. 8) used is 
dU: 
a, > EE (Z5 +t 2,1}. (IIL. 3) 


The vertical grid structure was introduced by Lorenz [1960] for use in 
the balance equations. 
All of the dependent variables occur at each x grid point. The x 


derivatives were approximated as 


(32), = Sie!) cre 


yn and (III. 4) 


al — (Uint Je Uj) Oe +S;) a (u; ust) (S; #S;.1) 


aX 4.Ax ° (III. 5) 


flere i Genotes the grid index and Ox the mesh length. This is a variant 
of the finite difference scheme that Arakawa [1966] developed [see also 
Ly, 1965]. The Arakawa scheme is designed to avoid the non-linear 
computational instability which was discovered by Phillips [1959]. In 


Prriyrcularstie nite difference approximation to the integral 
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Z race S3s9))\ ; (2), cage 


teezero when (III. 1) and (111.5) are used. The integral lis always 


é and @2 


zero in the continuous case. If the domain averages of uv’, V 
ame conserved by the non-linear terms throughout a numerical integra- 
tion, then no non-linear computational instability can occur. When 
finite differences in time are employed, these quantities will not be 
identically conserved and instability may arise. 

The continuous equations conserve the total energy E, which is 
given by Equation (II.21), If non-linear computational instability 
@eours, I cannot be conserved. It was found that EF remained nearly 
constant in the numerical solutions in this experiment, which indicates 
that the spurious computational production of energy was small. This 
was important because it was desired that the computational scheme 
would not introduce a smali-scale energy source which could produce 
a false jump. 


Centered time differences were used throughout this experiment 


with the exception of the first time step which was a forward time step. 


iS 








IV. ocCALE ANALYSIS 


A scale analysis is now performed that will indicate the dependence 
of the solutions on the Rossby number and Froude number. The inde- 
pendent variables are written 

+= bel se ie ane eee - (IV. 1) 
The dependent variables are scaled as follows: 

ne ee = Vy = Wwe = UIOINGe ea ad 


6 


(UNO,/g)0'. (IV. 2) 


The Rossby number and Froude number are defined as 


Ro = U/f{L and F = U/DN (IV. 3) 
mm@eto Neis @etined 26 
= ., | 
2 9@ {2 
Se cle 
n=[2 o (IV. 4) 


The time scale T and the v-scale V are functions of Ro and F while 
Pe w-scale is W = UD/L. 

An examination of the case where F£1 and FS Ro is now given. 
This should parallel the results obtained by Houghton [1969]. The 
appropriate time scale is T = L/DN, and the v-scale is V = (F/Ro)U. 

dV 


The later was obtained following Houghton by balancing 3 with fu 


in Equation (II.6). The non-dimensional equations for this case are 








Su y(u'u’) (uw) »¢’ eae (IV. 5) 
pe aes : ol 4 ie k — 








dv’ d(ulv’ vw ' (TN ao} 
Je 5 5 | dy! mi ya aE U = 0 


0 (48) (w'8) (IV. 7) 
20 ye | Oy SO + wis, 

eee ety a (IV. 8) 
dx’ dZ° 

ae Lg’ e 0. (IV. 9) 
oz! 


Hydraulic jumps are formed through the action of the non-linear terms 


7 at eae ef [eS Tee’ = a ak Vgc 
Das duypuceta lio Oe Perera LV a: ‘fa 


Ewen when these terms ait 
small (Fé 1), they can produce a jump ina non-rotating system. If 


a hydraulic jump is to be prevented, the Coriolis term in (IV.5) must 


dominate the non-linear advective terms, which implies that 


ee | IV. 10) 
Ri? E , ( 


The curve dividing the jump region from the no-jump region should 


then have the form 


F = AR: Gy eb 


where the constant A should be appreciably Jarger than 1. Houghton's 
analytical curve corresponds to A™ 6.5; Williams and Hori [1970] 


found that A = 6.0 to 7.5 was in agreement with their numerical 


\7 





solutions. As in Williams' and Hori's experiment, with A in the range 
6.0 to 7.5, the above scale analysis does not apply for F~w 1, since F 
is then greater than Ro. 

Iwo ©esions in| Ro are of interest whens» 475) elt ore lee 
terms in Equations (1V.5), (IV.6) and (IV.7) are of the same order 
and a jump would be expected to form. However, when Ro 1, the 
appropriate time scale is I = f-! and the v-scale is V = U. The non- 


dimensional equations then become 


- 








+ Ro b= zr re + a ie 0, (EV es 
dV" 5 (u'v’) d(v'w’) : 


plane with Equations (1V.6) and (1V.9). If the terms in Ro are neglected, 
the resulting equations describe an inertial oscillation in u and v; the 
Bemoperature field is independent of time: 

The scale analysis indicates that, when F is small, the jump 
boundary is given by F = tae and that as F approaches l, the 
boundary will be between Ro =1.0and Ro=0.1. These results were 


rested by direct numerical integration of the basic equations. 
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Vwi CONDETIONS 


The initial x-averaged stability is independent of z. A wave 


oa” 
nlc 


sinusoidal in x, the amplitude of which is directly dependent on the 
Froude number, is superimposed upon this basic state. 


The initial fields for this experiment are given by the equations 


u(x,0) = U cos(kx) i cos E (Z - 2) - 2 


v(x,0)= 0 and Sale 


Q(x, 0) = ae 


The constants which define the initial state have the following 





cos (kx) sin] 2 (2-2) -£2 (2-8), 


numerical values: 


D = 3km., 
a2 4.0 °K km. and (V.2) 


¢ 


Z - 0,03269 m see" OK 


The wavelength of the disturbance, L, is equal to 1,000 km. 
miter finite diiferences are introduced, it is necessary to adjust 


the values of u which are given by (V.1) in order that 


bez du; 
4S 5 see. a Vv. 
2B Bag -Fay)e0 v3 


PSI— 


This merely states that the finite difference vertical velocity at the 


upper boundary must be zero. 








The following initial conditions were also examined: 


u (x, 0) U cos(kx)cos (3) , 
v(x,0) =O and 


UNO, 
6(x,0) “Frc cos (kx) sin (2) + $5 (2-2). 


tt 


(V.4) 


} 





When f = 0 they lead to an exact finite amplitude wave solution 


which moves with the linear gravity wave speed. In this solution the 


non-linear terms in Equations (II. 7) and (II. 18) cancel out identically 


In this case no jumps can occur although they were observed for f # 0. 
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VI. NUMERICAL SOLUTIONS 


The numerical solutions are presented in non-dimensional form. 
The time for all results presented is given in units of the time scale, 
L/DN. Most of the results are for Fw 1. Enough results are pre- 
memted for F < ! in order to determine the jump curve. This will 
provide a basis for comparison with the work done previously by 
Houghton [1969] and Williams and Hori [1970] in the one-layer, homo- 
peneous system. 

eeduantity useful inthe andlysis of the @unap criterion 1s dietined 


( Uwe r Unsin) 


(\ime by 


The quantity ae is approximated by (u;4] - u-)/ Ax. Henemhe oui 
scripts max and min refer to the maximum and minimum overall 
points ata giventime. The non-dimensional quantity d is a measure 
of the width of the zone through which u varies most rapidly. It will 
be an overestimate of this width when the values of u at the edges of 
ile Zone are not near Mey ceed ee However, if a discontinuity 
Ewch as a pressure jump does occur, d must approach zero, since it 
Is inversely proportional to the gradient of u. 

Figure 2 shows the variation of d as a function of time (time unit 
iy ON) for F = 1.0 and Ro =a. Results are shown for three values of 


Ax = 27 L/n corresponding to n= 20, 40 and 80. As can be seen, d 


decreases to a minimum at approximately DNt/L = 0.4. 


Z| 








sh 7=X JO ONTeA JUSTSTI_ e& 0} spuodsaiiod aAano yea a 
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Detailed solutions are presented for three cases in which F = lI. 
In these numerical solutions Ax =2Z727tU/n, wheren= 80. Figure 3 
shows the time evolution of the non-dimensional variables u' = u/U 
mae y/U and 6" =1(0 = 6) /(UN6,/g) [where ® represents the horizontal 
average| for Ro=1. The u and 6 fields show that a hydraulic jump is 
forming at DNt/L = 0.37. The v field, which begins at zero, grows in 
opposite phase to the u and 6 fields through the action of the Coriolis 
fence. Once the v field is present, the Coriolis force will tend to 
reduce the amplitude of the u field while the amplitude of the @ field 
remains relatively unchanged throughout the period of integration. For 
this value of the Rossby number the Coriolis force was unable to 
prevent the formation of the hydraulic jump. 

Figure 4 contains the evolution of the non-dimensional variables 
u’, v' and 9' for Ro =0.1. It appears that an inertial oscillation is 
present with the 6 field remaining very nearly constant throughout the 
integration period. It was found that when the integration time was 
eesended, though, the fields never were restored to the original as 
was shown for the one-layer, homogeneous case by Williams and 
Hori [1970]. 

The results for Ro = 0.35 are presented in Figure 5 using the 
feiaie format as above. In this case there is no evidence of any 


meauction in tne u field as the jump time is approached. 


es 
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Results for F = 0.04 and Ro = 0.09 are contained in Figure 6 
showing the evolution of u', v' and @' for a typical jump case when 
F€1. Figure 7 contains the evolution of the variables for F = 0.04 
and Ro = 0.07 giving the results for a non-jump case when F is again 
much less than l. 

In Figure 8 and Figure 9 the results are presented for F = 0.3 
oma Ro = 0.3 and 0.1, respectively, showing the time evolution of 
u', v' and @' when F is in the range 0.1 to 1.0. 

The above figures are taken from the results of the lowest layer 
in the model. Figure 10 is included in order to show what is taking 
place in the other layers at time zero and at an approximate jump 
mac. The case presented is for F = 1.0 and Ro=1.0. As can be 
seen, only the lower half of the model is shown. This provides a 
picture of what is taking place in the most active of the layers in the 
maodael. 

ine d procedure was Used to determine the curve separa. 
jump cases from the non-jump cases over the range F = 0. 04 up to 
We 0. These results are shown in Figure 11, with the logarithm 


2 and 


of F plotted against the logarithm of Ro. The lines F = 6.0Ro 
F= 7. sree Are Included ton recOnvenience, =simese Teculis tore quation 
(IV. 11), which was derived by scale analysis leading to the conclusion 


miat A is in the range 6.0 to 7.5. This relation does not hold for Pa i, 


but this was to be expected since the analysis required that FS Ro. 
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The non-jump solutions included in Figure 1] are determined over 


an integration time of OJ SBINE 


ot 








VII. CONCLUSIONS AND RECOMMENDATIONS 


The formation of hydraulic jumps ina rotating, continuously 
stratified fluid has been investigated. A scale analysis indicated that 
jump cases should be separated from non-jump cases by the curve 
F = ARo“ in the F-Ro plane. (ihe Nossmyenunim .e) i Olamdetnc wer ouc 
number F were utilized as direct inputs to the equations and used as 
a means of determining the initial scales (U = FDN). It was determined 
that the results of this experiment (with F appropriately redefined) 
closely paralleled the results of Houghton [1969] and Williams and 
Hori [1970] in which one-layer, homogeneous fluids were studied. 


e Pie J e 


ro ce + % a=— = Pe lm ae Nees ‘~os . x ~ ie 2°" a4 TAN 1 = ? 
wee Trasiitc Matar SULUciaGgme Cl Using SaaS 12a F mplv far; ofS 1 Rane 


range 6.0-7.5 and that the relation F = IN Seeye 


holds up to at least 
Fe 0.35. Extrapolation to F = 1.0 showed jumps occurrimge for 
smaller values of Ro than were given by the above curve. ‘Thus it can 
be seen that the jump criterion becomes more independent of Fas 
F = 1.0 is approached. 

This work should be extended to an atmosphere in which the jump 
region is bounded by a stable layer rather thana rigid plane as in 
this case. Sucha study could require additional layers with further 


consideration given to layers of varying thickness. Further considera- 


tion should be given to moisture in future experiments. 
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